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Abstract: 

Solid state quantum condensates often differ from previous examples of condensates (such as 
Helium, ultra-cold atomic gases, and superconductors) in that the quasiparticles condensing have 
relatively short lifetimes, and so as for lasers, external pumping is required to maintain a steady 
state. On the other hand, compared to lasers, the quasiparticles are generally more strongly 
interacting, and therefore better able to thermalise. This leads to questions of how to describe 
such non-equilibrium condensates, and their relation to equilibrium condensates and lasers. This 
chapter discusses in detail how the non-equilibrium Green's function approach can be applied 
to the description of such a non-equilibrium condensate, in particular, a system of microcavity 
polaritons, driven out of equilibrium by coupling to multiple baths. By considering the steady 
states, and fluctuations about them, it is possible to provide a description that relates both to 
equilibrium condensation and to lasing, while at the same time, making clear the differences 
from simple lasers. 
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Bose-Einstein condensation (BEC), the equilibrium phase transition of weakly interacting 
bosons, was realised over ten years ago in ultra-cold atomic gases. After long and strenuous 
efforts to observe this state in solids, BEC of polaritonsf]]] and of magnonsfil^ were reported. 
These reports followed observations of related effects for excitons in quantum Hall bilayers[3|, 
spin triplet states in magnetic insulators 0| and excitons in coupled quantum wells [5- 7]. Solid- 
state condensates depart from the archetypal BEC in several ways. Most importantly they 
live for short times and rely on external pumping. Indeed, it was the decay, and consequent 
lack of equilibrium, which for a long time presented the obstacle the realisation of solid-state 
BEC. Even if one can accelerate thermalisation, the decay, and the consequent flux of particles, 
remains a more important effect in solid state than it generally does in cold atomic gases, or in 
other quantum condensates such as superfluid Helium. 

When considering whether such a system may be treated as equilibrium or not, there are 
several distinct characterisations of the degree to which the system is non-equilibrium. The most 
obvious compares particle lifetime to the time required for collisions to thermalise the system, 
determining the extent to which a thermal distribution may arise. The timescale for establishing 
a thermal distribution within one part of the system can however be quite different to that 
for establishing either thermal or chemical equilibrium between different parts of the system. 
Another characterisation of whether non-equilibrium physics is relevant arises from comparing 
the linewidth due to finite particle lifetime to the temperature of the system, thus determining 
whether lifetime or temperature effects dominate coherence properties. Table [TrT] summarises the 
typical timescales and energy scales connected with different examples of metastable quantum 
condensates. It is clear that the ratio of thermalisation time to the particle lifetime is generally 
somewhat larger for solid-state condensates than it is for cold atomic gases. If one instead 
compares the ratio of the linewidth due to decay to the characteristic temperature, polaritons 
stand out as having a decay linewidth of the same order of magnitude as their temperature. 
As such, polaritons are good systems in which to study effects of finite lifetime on coherence 
properties. 





Lifetime 


Thermalisation 


Linewidth 


Temperature 


Atoms [8] 


10s 


10ms 


2.5 X 10-^^meV 


10~**K lO-^meV 


Excitons [9] 


50ns 


0.2ns 


5 X 10~^meV 


IK O.lmeV 


Polaritons f 10] 


5ps 


0.5ps 


0.5meV 


20K 2meV 


Magnons[2j 


1//S 


100ns 


2.5 X lO-^meV 


300K 30meV 



Table 0.1: Characteristic timescales and energies for: particle lifetimes, times to 
establish a thermal distribution, linewidth due to finite lifetime, and characteristic 
temperatures for various candidate condensates. Comparison of the first two de- 
scribes how thermal the distribution will be; comparison of the later two determine 
the effect of finite lifetime on coherence properties. 



Because, as we will discuss further below, polariton condensates provide such a clear illus- 
tration of the properties of non-equilibrium condensation, we will focus on them in particular. 
Microcavity polaritons are the quasiparticles which result from strong coupling between pho- 
tons confined in a semiconductor microcavity, and excitons in a quantum well. By changing the 
detuning between the excitons and photons, and by changing the strength of an external pump 
that injects polaritons, one can modify the polariton mass, density and the effect of interactions 
between polaritons. A more detailed introduction to microcavity polaritons and semiconductor 
microcavities can be found in several review articles and books Ill4l7] 

The intrinsic non-equilibrium and dissipative nature of solid-state condensates, especially 
of polaritons, brings connections to other systems exhibiting macroscopic coherence, i.e lasers. 
With the realisation of more complex, interaction dominated lasers, such as random lasers (see 
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e.g. Refs. H, 3]) or atom lasers (e.g. Refs. 13- ll]), this connection is particularly pronounced. 
Compared to simple lasers, polaritons are however more strongly interacting, and therefore much 
better able to thermalise than are photons, and so in many ways solid-state condensates can be 
viewed as being somewhere in between an equilibrium BEC and a laser. At the same time, at 
large temperatures and/or in the presence of large decoherence mechanisms and large pump- 
ing the same microcavity system supports a simple lasing action. In this context, microcavity 
polaritons provide particularly excellent playground for studying coherence in a dissipative en- 
vironment, and the differences and similarities between condensates and lasers. Clearly, an 
approach which takes into account the non-equilibrium and dissipative nature of this new state 
of matter, as well as strong interactions, multimode-structure, low dimensionality and finite size 
is necessary. 

This chapter will discuss a theoretical approach to modelling quantum condensates that are 
driven out of equilibrium by a flow of particles through the system. We therefore consider cou- 
pling the system to baths, which can transfer energy as well as particles to and from the system. 
With such baths, we find that the behaviour of a simple laser can be recovered in the limit of high 
temperature baths. A different scenario of how decoherence affects condensation can be found if 
one considers static disorder — i.e. allowing scattering, but with no transfer of energy to or from 
the system. Such a problem 23j, |2j] is closely related to the Abrikosov-Gorkov approach to dis- 
ordered superconductors 25|. As in the case of superconductors, one finds a distinction between 
"pair-breaking" and "non-pair-breaking" disorder (respectively magnetic and non-magnetic im- 
purities in the superconducting case). As expected from Anderson's theorem [26|. the coherence 
associated with the condensate leads to a gap in the exciton density of states, which makes 
the condensate robust to non-pair-breaking disorder. With pair-breaking disorder, decoherence 
eventually destroys the gap and finally the condensate, but for small amounts of decoherence, 
the gap protects the condensate. A similar scenario also exists in the ultra high density limit, 
where excitons are destroyed by screening, leading to an electron- hole plasma phase[27[, which 
can nonetheless support lasing. While we focus in this chapter instead on the effects of particle 
flux, and baths that can transfer energy, these other results illustrate that there are a variety 
of ways in which decoherence can either suppress or modify the properties of a condensate. In 
principle one can have both a crossover from a polariton condensate to a regular laser (weak 
coupling but still excitonic gain medium, as discussed here), and a crossover to a particle-hole 
laser (weak coupling, electron- hole plasma, if screening is strong). 

The approach to modelling the condensate with a flux of particles presented in this chapter is 
based on work by the authors in Refs. [2S. ,29]. While in those works, the results were derived and 
presented making use of the non-equilibrium path integral approach 30l|. both the results and 
their theoretical basis can be understood without this technical background, by considering the 



diagrammatic approach to calculating non-equilibrium Green's functions [3ll433t| . The particular 



aim of this chapter is therefore to review some of these results, illustrating in some detail how 
a steady state non-equilibrium system which develops spontaneous coherence, can be treated in 
the non-equilibrium diagrammatic formalism. At the same time, this approach will provide a 
natural language to highlight the way this system relates both to equilibrium condensates and 
to lasers, and to understand the ingredients that makes it differ from these limits. 

There are a number of other known approaches to describing systems driven out of equi- 
librium by coupling to multiple baths. Those that have been applied to microcavity polari- 
tons include: quantum kinetic equations [34 - 41 1 . Heisenberg-Langevin equations [43], stochastic 
methods for density matrix evolution (i.e. truncated Wigner approximation) [43. ,44] : as well as 
mean-field approaches, considering the complex Gross-Pitaevksii equation, in some cases includ- 
ing also coupling to reservoirs or thermal baths 4^47|. While this chapter does not intend to 
review the merits of each of these approaches, it is worth noting that in general, these approaches 
are all connected. The connections between many of them can simply be seen by looking at 
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their relation to the non-equihbrium diagrammatic approach. As discussed in [33l . |48| . the 
quantum Boltzmann equation can be derived as an equation for the distribution function that 
appears in the Keldysh Green's function, along with a Wigner transformation from F(t,r,t' ,r') 
to F{T, R,uj,p). It will become clear from the discussion in Sec. El that there is a close analogy 
between the Keldysh Green's functions and the Heisenberg-Langevin equations, with the bath 
Green's functions describing the same physics as the correlation functions of the bath noise 
operators in the Heisenberg-Langevin approach. There also exists a connection between the 
approach described here and density matrix evolution. The single particle density matrix is 
given by {ip^ {r,t)ip{r' ,t)), and thus corresponds to an appropriate combination of equal time 
Green's functions. The density matrix naturally gives single time expectations of appropriate 
observables, it is also possible to derive two-time correlations from the time evolution of the den- 
sity matrix, by means of the quantum regression theorem [49]. The quantum regression theorem 
however relies on making an additional Markov approximation regarding the bath occupations, 
as well as a Markov approximation for the bath density of states (5o|. The Keldysh Green's 
function approach does not require this additional Markov approximation; and in fact Sec. H] 
will show how making this further approximation restricts the conditions for condensation to 
occur. 

In order to illustrate the application of the non-equilibrium technique, we consider a specific 
model of microcavity polaritons, starting from disorder localised excitons strongly coupled to 
cavity photons [5 1!-[53|. In this model, interactions between excitons are included by treating 
the excitons as hard-core bosons, allowing one exciton, but no more, to occupy a given disorder 
localised state. For the discussion presented here, using this model provides a number of technical 
advantages: it connects closely to the idea of gain from two-level systems that is typically 
used in models of simple lasersfi^, making the comparison to lasing straightforward; and it 
automatically includes nonlinearity of the excitons, allowing this nonlinearity to be described by 
the properties of the exciton representation, rather than requiring higher order diagrammatic 
corrections. In addition, in an equilibrium situation, the mean-field theory of this model is 
known to give a reasonable description of the critical temperature, except at very low densities 
where fluctuation corrections become importantjs^. 

This chapter is organised as follows; section [1] introduces the model Hamiltonian, and its 
coupling to baths. Section [2] then describes the approach we will take to modelling this system, 
reviewing some standard results of the non-equilibrium diagrammatic technique that will be 
used later, and discussing the mean-field approach we use to find the steady state. In order 
to evaluate this mean-field condition, it is necessary to determine the effects of the baths on 
the system, by calculating particular self energy diagrams, these self energies are presented 
in section [3l Section U] then discusses the mean-field theory, considering how it can recover 
both equilibrium results in one limit, as well as the description of a simple laser in another 
limit. Section [5] discusses fluctuations about the steady state, analysing stability, and further 
illuminating the connection to (and distinctions from) a simple laser; section then provides a 
more qualitative discussion of the fluctuations of the condensed system, focusing in particular 
on the combined effect of finite size and finite lifetimes. 



1 Polariton system Hamiltonian, and coupling to baths 

As explained above, we consider a model of excitons as hard core bosons coupled to propagating 
photons. To write the Hamiltonian for hard-core bosons, it is convenient to introduce fermionic 
operators i| , a| , such that the two fermionic states represent the presence or absence of an 
exciton on a given site, hence the operator bla- is the exciton creation operator. With this 
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notation, the system is described by: 

i^sys = Yl " "1"^) + E^^^^^fe + T^^M^^ + H.C.), (1) 

i k i,k 

where ej is the energy of a locahsed exciton state, gi is the exciton-photon couphng strength, 
and Wfc = Wo + /c^/2mphot is the dispersion of cavity photons. As sketched in Fig. 10.11 this 
can then be driven out of equihbrium by couphng to two baths, so the the system evolves 
under the full Hamiltonian H = Hgys + -ffsys,bath + -f^bath- It will be useful later on to divide 
the coupling to baths into coupling to the pumping bath, and coupling to the decay bath 
-ffsys,bath = -f^sy^s^ath + "^syT.bath where the forms of the coupling to the pumping and decay baths 
are: 

KZlt^ = E («! A„ + blB.^ + H.c.) , Fj-^,, = ^ C,,k + H.C.) . (2) 

n,i p,k 

Here Tn,i is the coupling to a pumping bath, described by the fermionic operators Bn,An, and 
Cp,k is the coupling to decay bath, describing bulk photon modes ^'p. The bath Hamiltonian is 
taken to have the simple quadratic form: 

^bath = E {^nK - AiA^) + E ^P^P^P 
n p 







System 




Bulk photon modes 


Pumping Bath 








— mode 

















Figure 0.1: Cartoon of system, consisting of photons strongly coupled to excitons, 
and external baths, describing pumping and decay. Adapted from Ref . (29| . 



Describing the pumping reservoir and photon decay as baths means that we assume these 
both contain many modes (i.e. are much larger than the system), and thermalise rapidly com- 
pared to the interaction with the system. These assumptions mean that one may impose a 
particular distribution function on the bath modes, and then determine what distribution the 
system adopts; we will take a thermal distribution for the pumping bath, specified by a bath 
temperature and chemical potential, and we will assume the bulk photon modes are unoccupied. 
Note that we do not explicitly introduce any system chemical potential, as the density of the 
system will be fixed by the balance of pumping and decay, however a natural definition of the 
system chemical potential will arise later. 



2 Modelling the non-equilibrium system 

The Keldysh non-equilibrium diagrammatic technique[31] is an approach well suited to dealing 
with the kind of non-equilibrium steady state which we consider here. Section [5] briefly sum- 
marises the concepts that will be important in the remainder of this chapter; for a more complete 



introduction, see for example Refs. [3ll433l|. Within this diagrammatic approach, we will then 
determine the possible steady states of the system by a mean-field approach, introduced in 
section [2j 



2. MODELLING THE NON-EQUILIBRIUM SYSTEM 



5 



Non-equilibrium diagram approach 

In order to determine both the spectrum (i.e. the ground and excited states, taking into account 
interactions and couphng to baths), and the non-equihbrium occupation of this spectrum, it is 
necessary to calculate two linearly independent Green's functions; it is convenient to make these 
the retarded and Keldysh Green's functions: 

D^it,r) = -ie{t) (^m,r),^P^O,0)]^) , (t^r) = -i (^m,r),^p\0,0)U) . (4) 

Here, [ip, ip^^ indicates the commutator (anti-commutator) of ^ and ip^ . These Green's functions 
can be written as time-ordered products of fields by introducing the Keldysh contour, shown in 
Fig. 10.21 Each point on this contour is labelled by {t, {/, 6}), where the /, b label whether it is 



b 



f 

Figure 0.2: Keldysh closed-time-path contour, which can generate multiple order- 
ings of fields 



on the forward or backward branch. We then introduce the contour time ordering Tc, such that 
fields on the backward contour are always later than those on the forward contour, and that 
pairs of fields on the backward contour should appear in reverse order. By then introducing 
symmetric and anti-symmetric combinations of these fields ■i/'i = [i^{tj f) ^ i^itjb)] /\/^, one 
may write the Green's function: 



D 







Tr 



ip-{t,r) 



(5) 



Here, refers to the advanced Green's function, which is the Hermitian conjugate of the 
retarded Green's function. 



Given the above time-ordered products, one may use standard methods SJ, l55| to write a 
diagrammatic expansion, by writing the Heisenberg picture fields in terms of the interaction 
picture fields ip{t): 

m = U-\t)mUit), m = e*^°Ve-*^°*, (6) 



where H = Hq + -ffint, and Hq is "free", meaning that it is simple to write expectations of 
products of fields evolving according to Hq. By formally solving the equation for U{t), one may 
then write the Green's functions in the following form: 



D 

U 



i(Tr 



|:<;;:5)(*i(o.o).^uo.o,)t, 



exp 



-i / Hi^t{t)dt 
'C 



exp 



Hintit,f)-Hintit,b))dt 



(7) 
(8) 



The diagrammatic expansion then follows by expanding the exponential, which produces ver- 
tices coupling free fields, and connecting these vertices by lines representing the Green's functions 
of the free fields. Compared to other diagrammatic expansions, the only extra complication is to 
keep track of the it labels on the fields, both in the matrix structure of Keldysh/retarded/advanced 
Green's functions, and in the form of U. 
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In the following, we will frequently make use of the Dyson equation [33|, |5J, |55|], D 



1-1 



-1 







S, and so it is useful to record the free inverse Green's function. The inverse Green's 



function has the structure: 



D 



-1 







-1 




K 



(9) 



where [D = — [D^] ^ [-^^] ^ - Using the results for a free field, one has 

[dI^]~^ =LO-Uk + iv, [D^'f = {2i7j)i2nB{u;) + l), (10) 

where t] is infinitesimal. All of the results noted above assume bosonic fields; the results for 
fermionic fields are similar, but commutators and anti-commutators are interchanged in the 
definitions of Keldysh and retarded Green's functions. 

For our particular model of microcavity polaritons, the division of the full Hamiltonian into 
Hq and Hint will be to take: 

Ho = Y, ^^ibl^ - + Yl ^-fc^I^fc + E 5^^o {e'^^'ajb^ + e-^^^'bla^) + i^bath (H) 

i k i 

where -00 is a mean-field coherent photon field, as discussed in the next section. This means that 
Hint will contain the system-bath interactions, as well as the interaction between the two-level 
systems and incoherent photon fluctuations. In the following we will however generally focus on 
one part of Hint at a time. 



Mean-field condition for coherent state 



For a system coupled to multiple baths, the mean-field theory can no longer be thought of as 
minimising free energy, but rather as a stable self consistent steady state. For a condensed 
solution, one looks for a steady state of the form (tpk) = '4'o^w{~'>'fJ'St)^k,o = 'ipo{'t)6kfl, where 
Us is introduced here merely as part of the steady state ansatz, but it will be seen to play a role 
analogous to the equilibrium chemical potential. To be a self-consistent solution, this ansatz 
must satisfy the Heisenberg equation: (idtip) = {[ipjH]), and so: 



(12) 



The expression {al{t)b-{t)) describes the polarisation of the two-level systems, and can be written 
in terms of the Keldysh Green's function, as: 



(13) 



As well as this self-consistency condition to determine the coherent field amplitude and the effec- 
tive system chemical potential the mean- field approach can also be used to give an estimate 
of the polariton density. This density will be used later in producing the phase diagram of the 
polariton condensate. The mean- field estimate of the total density is given by the combination 
of the photon density l^/'oPj and the fermion density (i/2)Tr[G/f — ]. 



3. EFFECTS OF BATHS ON SYSTEM CORRELATION FUNCTIONS 
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3 Effects of baths on system correlation functions 

In the above, we found that the mean-field condition could be written in terms of the two-level 
system Green's function, and the expectation of the decay bath fields. In this section we will 
discuss in detail the treatment of the baths and their effect on system's correlation functions, 
which will then determine the conditions under which a condensed solution may exist. Most 
of the effort, in Sec. El will be dedicated to finding {u) including the effects of pumping. 

Before doing this, section [3] will address the simpler problem of how {^p{t)) can be related to 
the decay bath Green's function and thus evaluated. 

Decay bath and (^p) 

To calculate (^p(t)) in terms of non-equilibrium Green's functions, one may first use the inter- 
action picture, in terms of the system-bath coupling, to write ^p(t) = U~^{t)4!p{t)U{t). Here, 
U{t) is the time-ordered exponential as in Eq. Q: 



U{t) = Texp 
Then, consider inserting a factor: 
1 = T exp 



sySjbath*' 



,,/fVdccay 
"■'^ -"sys,bathl'' ) 



■ Texp 









(14) 



(15) 



either before or after ^p{t). The resulting expression implies that one has: 

i^pit)) = {Tc[^pit,f)U]) = {Tc[^p{t,b)U]) = i=(rc[§p,+(t)C/]), 



(16) 



where the last equality has made use of the fact that if the expectation of ^{t,f) and ^(t, 6) 
match, then the expectation of ^-{t) must vanish. We are interested in particular in the value 
of this expectation (^p) when we consider the system in the mean-field approximation. In this 
case the system bath interaction term is given by: 



dtH' 



c 



decay 
sys,bath 



it) 



dt Cp,o V2 (0 V'o (t) + ro mp,- (t)] . 



(17) 



With vertices given by this interaction, the set of diagrams involved in evaluating Eq. (I16p 
is particularly simple: the only possible connected diagram is one with a single bath Green's 
function connecting the source term in -ffsys,bath to the field ^+ that we want to measure. As 
such, the sum appearing in Eq. (fT2l) can be written as: 



EC? 



p,0 



dt'D^, it,t')Mt') 



(18) 



The simple form this equation takes is also the form one would find by making the Born ap- 
proximation; i.e. assuming that Cp,o is small, so that terms like ^pCpfl should be kept, and 
neglecting terms involving any higher power of Cp,o- However, in the current case, because of 
the linearity of the coupling, no other connected diagrams exist, so no assumption of smallness 
is required in order to neglect higher order terms. 

For a free bath, one may write ^p{t) = e~*'^p*^'p, and so the bath Green's function is given 
by D^^ (t, t') = —i6{t — i')e~*'^p(*~* ). Taking a Markovian approximation for the bath density 
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of states and coupling [i.e assuming that the product of the bath density of states N^{uj) and 
the square of the system-bath coupling (p^ are constant: tt^p^qN'^ (uj) = k] then gives: 

Yl Cp.oe"'""^^*"*'^ = 2Ac<^(t - t'). (19) 
p 

Putting this all together, the net effect of the decay bath on the self-consistency equation is just 
to add a decay term, so the result may be written as: 



Pumping bath and 

The remaining task is to find the matrix of fermionic Green's functions in the four by four space 
resulting from the a, b fermionic fields and the it space associated with the closed-time-path 
contour. As above, we will take the interaction Hamiltonian to be the coupling between the 
system and the bath. This leaves the free fermion Hamiltonian: 

^0^"^' = E ^^(^1^^ - + E 9^^^ {^'''''^Ih + e-^-^'bla,) . (21) 

i i 

It is possible to diagonalise this Hamiltonian by a unitary transformation, and thereby write the 
appropriate free Green's function, however it is first necessary to remove the time dependence 
introduced by the form of the ansatz for the coherent field. This can be achieved by a gauge 
transformation: 



2 



E m - + E i^nBn - AiA, 



(22) 



such that b — )• 6e~*^*'*/2,a — )■ ae^^^^^'^, which removes the time dependence of the mean-field 
photon to fermion coupling. The gauge transformation for the bath modes that also appears in 
Eq. (j22|) is necessary to ensure no time dependence is introduced into the system-bath coupling 
terms. The net result is to replace — )■ = — //5/2 in Hq^^, and to shift the bath Green's 
function in frequency by zizns/2. 

After the above transformation, the Hamiltonian can be diagonalised by the unitary trans- 
formation: 

bi \ _ f cos{9i) sm{6i) \ [ l^i 



ai J \ -sm{9i) cos{9i) J \ Oi J ' ^'^^^ 

after which the free Hamiltonian takes the form Hq^^ = Ei{(3jP- — a^a^), where tan(20j) = 
—giip/ii and Ef = e? + gfipQ- Since the Hamiltonian is diagonal in the /3,a basis, the retarded 
Green's functions in that basis are just [u ^ Ei + ir/]"^ (where r] is infinitesimal), and so the 
retarded Green's functions in the b, a basis can be written as: 

j^, _ / cos6'i sin^i \ / [u - Ei + ii]]^'^ \ / cos6'j -sin6'j \ 

~ V -sin^i cose, J \ + E, + ir]]-^ J \ sin Bi cosOi J 

{u + ir/)2 - Ef V 9ii^0 V -€i+iri ) 
\G^]-' = 1^-^^ + '^ y (25) 



3. EFFECTS OF BATHS ON SYSTEM CORRELATION FUNCTIONS 
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Just as for the free bosonic Green's functions described in Eq. (fTO|) . [Gq"^]^^ is infinitesimal. 
Since coupling to the pumping baths, which is discussed next, will add a non-infinitesimal 
Keldysh self energy, we will neglect this infinitesimal contribution. Therefore, the expression for 
[Gq~\ and its Hermitian conjugate are all that is needed of the free Green's function. What 
remains is to determine the self energy that arises from coupling to the pumping bath 

Taking the part of the interaction Hamiltonian due to -f^^™ath' ^^'^ inserting it into the 
definition of U in Eq. ([8]), one has that the interaction vertices are generated by: 



c 



"''-"sys.bath 



dt 



— oo 
oo 



dt 



dlit,f)A^{tJ)-dlit,b)A^it,b) + 
«LW4- W + alW^n+W + • • • 



(26) 



As each pumping bath couples either to only a modes or to only b modes, there are no off 
diagonal self energy terms in the a,b basis. One may therefore concentrate first on S^t^, and 
Sfjtft will follow by analogy. Written as a matrix in the Keldysh space as defined in Eq. ([7]), one 
has: 

y++ y+- \ 

^ (27) 



at a at a 
aTa aTa 



The lb labels determine the label of the incoming/outgoing fields, and it is clear from Eq. ()26p 
that + fields couple to — bath fields and vice versa. Thus, an example self energy diagram is: 

giving the equations 



a 


A a 




n 




= ^ ^l,nG 

n 


SV(i,t') 


= X] ^InG 

n 


s--(t,t') 


= X] ^InG 

n 



At A 



At A 



At A 

++ 
At A 







n 

+i ^lJ{t' - t)e+*'^" = +i-i5{t - t') 

n 



In the last three lines, the Markovian limit has been taken to give the final equality. In the last 
line we have used: 



FA{t) 



dv 
— f 
27r 



-iut 



FA{y) 



FA{v) = l-2nA{v) 



(28) 



where the form of the distribution function F comes from the form of the equal-time Keldysh 
Green's function F^(f^) = {AnAn — AnAn). As a function of frequency, the self energy matrix 
in Keldysh space is thus: 

^7 
-i7 —2ijFA{i^) 



(29) 



The matrix for S(,tb(i^) is identical except that Fa{i^) — )• Fsiv)- 

Combining the free Green's function and self energy, we may write the entire inverse Green's 
function in the basis (6+, a+, 6_, a_) as: 



( V -ii-i-i -giipo \ 

-gii'o v + li- 

u-ii + i-f -Qiipo 2i-fFB{i') 

-Qii^o ly + ii + n 2i-^FA{v) J 



(30) 
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Clearly, if is zero then the a,b fields decouple as expected. However for non-zero "00! the 
Keldysh Green's functions of the two fields get mixed, so that their occupation is set by a 
balance of the pumping and the effects of the coherent photon field. 

To complete our analysis, we should invert the above matrix to find the Keldysh block, and 
the a^b component of that block. Using the Keldysh block structure of Eq. ([9]), and in particular 
that = -G^iG-^G^, one may write: 



G^M - 



[{u + ijy - Ef][{u - ij)^ - Ef] 
-ii + ij -gitpo \ [ Fb{v) \(v-ei-i-f -Qiipo 



giipo u + ii + i-f J \ Fa{i^) J \ -giipo v + ii - i^i 



(31) 



For the mean- field condition in Eq. ([El), we require in particular the component which 

has the form: 

' [(.'-Ei)2 + 72][(i. + E,:)2 + 72] ^ ' 



4 Mean-field theory and its limits 



Putting the component of Eq. (|32l) into Eq. (I20p gives the self-consistency condition (equa- 
tion for the condensate) of the mean-field theory: 



/ .^, V- 2, / du [Fb + Ea)v + [Fb - FA)[ei + Z7 

{u:o -l^s- = 2^ ft ^07 / ^ ^ _ ^ ^ ^ • (33) 



2, f dv {Fb + FA)y + {Fb - FaWi + i7) 

Equation ()33p is central to our analysis. This equation is rather powerful, in that it combines 
several well known theoretical results within a single framework. As will be shown in section [H 
in the equilibrium limit (where the system-bath couplings are taken to zero) Eq. (|33l) reduces 
to the gap equation which applies throughout the BCS-BEC crossover. In the opposite highly 
non-equilibrium limit (see section H]) it reduces to the standard laser condition. At low densities 
it reduces to the (complex) Gross-Pitaevskii equation, discussed in section HI As such, this ap- 
proach highlights the connections between these apparently different descriptions of condensates 
or lasers. 

The functions Fa,b appearing in Eq. were defined as Fa,b = 1 — 2n^^B, where nA,B 
are bath occupation functions. These occupations are taken to be externally imposed, and can 
be chosen to have any form relevant to a particular physical situation. Here, we will choose 
these to be thermal and at equal temperatures but different chemical potentials. Noting that 
the fermionic states were supposed to represent two-level systems (or excitons), we take the 
occupations to satisfy ua + nB = 1. This therefore requires that we have: 



Fa,b{i^) = tanh 



1^ ( .IJ'B-IJ'S 
— [ V ± 
2 



(34) 



where /i^ is an adjustable pumping bath chemical potential, and appears in this expression 
due to the shift arising from the gauge transformation in Eq. (|22l) . Schematically, this situation 
is illustrated in Fig. l0.3| one can see that F^(— e)+Fe(e) = 2[1— nA(— e)— ^^(e)] = 0. Physically, 
this pumping process is most closely related to electrical pumping. Note that, in the absence of 
any other processes, contact between the two-level systems (excitons) and the pumping reservoir 
would control the population of the two-level systems, and so one would have: 



(6^6 — a) a) = nB(e) — nA(— e) = —tanh 



P ( I^B , IJ-S 

2 l^-l- + T 
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Thus, by pumping with a thermahsed bath, one will find a thermalised distribution of excitons. 
Therefore, in the context of polaritons this pumping scheme resembles closely pumping from a 
thermalised excitonic reservoir, which is often the case in the experiments. 




Equilibrium limit of Mean-field theory 

The simplest limit to recover from the non-equilibrium self-consistency equation is that of ther- 
mal equilibrium. This corresponds to taking 7, k — )• 0. Since the self-consistency equation 
included only the coupling between mean-field photons and the decay bath, there is no way that 
a thermal distribution can be set by the decay bath. On the other hand, the pumping bath can 
set a thermal distribution, so to recover a non-trivial equilibrium distribution one should take 
K — 7- first, and then 7—7-0. If k = 0, then the imaginary part of the right hand side of Eq. ()33p 
must vanish. In order to satisfy this, without restricting the range of solutions of the real part, 
one must choose Fb{i^) = Fa{i^)- In terms of the distribution functions written in Eq. (I34p . 
this clearly means ns = MB- Physically, this means that in the absence of decay, the chemical 
potential of the condensate matches the pumping bath. 

After fixing /i5, the remaining part of the equation becomes: 

2 ; /" 2tanh (/3z^/2) 

We may then take the limit of small 7, by using: 

277^ 27r 



hence we find: 



(a;o - /xb)Vo = di^tanh [6{u - E,) - 6{u + E,)] 
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This is the equilibrium result j5ll453l]. but with the two-level constraint on the fermions imposed 
only on average. ^ Note, that this is the standard mean-field gap equation of the BCS-BEC 
crossover theory [sH ]. 

High temperature limit of Mean-field theory - simple laser 

The opposite extreme to the equilibrium condensate is the limit of a simple laser, which can also 
be recovered from Eq. ([33l) . Before showing how this limit can be recovered from our theory, 
we first provide a brief summary of the threshold condition of a simple laser, and express it in 
similar language to the above self-consistency condition. The equations describing the steady 
state of a laser can be derived starting from the well-known Maxwell-Bloch equations: 

dfipo = -ioJoV'o - KTpo y^ffjPj, (39) 

i 

dtPi = -2ieiP - X^Pi + gi^PoNi (40) 
dtNi = A||(iVo-A^i)-25i(^o^i + ^>o). (41) 

These equations can be understood as originating from considering a Hamiltonian like Eq. ([1]), 
with Pi = —i{a\hi) , Ni = {h\bi — a\ai) . One then writes the Heisenberg-Langevin equations, with 
a Markovian set of baths distinct for each two-level system, and then takes the semiclassical ap- 
proximation to drop bath noise operators. The value Nq is the bath inversion imposed by the 
pumping. Note that with coupling to such a Markovian pumping bath, there is a discontinuous 
jump between the allowed steady states with no decay, and the laser-like solutions found for 
any non-zero pumping|57|. In particular, with pumping and decay, inversion is always required 
for a condensed solution of these Maxwell-Bloch equations, so they cannot smoothly interpolate 
between a condensate and a laser. Such behaviour should not be too surprising, as a frequency 
independent (Markovian) bath occupation corresponds to an infinite temperature, and so even 
arbitrarily weak coupling of the system to an infinite temperature reservoir may destroy the 
condensate. With more realistic models of pumping, such a discontinuous jump need not nec- 
essarily occur. One should thus interpret the microscopic origin and consequent behaviour of 
Eqs. (|39l - HT|) with some caution. However, since Maxwell-Bloch equations of the above form 
are frequently used as a simple model of a laser, it is instructive to see what approximations 
they would correspond to in terms of our non-equilibrium formalism, in which the microscopic 
description of the pumping is better controlled. 

Starting from these Maxwell-Bloch equations, the self-consistency condition for a macro- 
scopic photon field i^oit) = -i/^oe"*^* can be written as: 

{-i^ + iuo + K)ipo = ^ giPi, {-ifi + 2iei + X±)Pi = giipoNi, (42) 

i 

which can be combined to write a single self-consistency condition: 

{ivo - n - iK)i>o = - y — — . (43) 

^—^ lei — iX±_ 



^Imposing the two- level constraint on average, the equilibrium expectation of the inversion (6^6 — a^a) can be 
written as: 



1 + e'S-B + e-/3^5 + 1 (e,8i5/2 _^ g-/3i;/2)2 ^^'''^ 2 J " ^^^^ 

Were the two-level constraint imposed exactly, the result would instead be: {e^^ — e^^^) j (eP^ -I- e~^^) — 
tanh iPE), as the zero and doubly occupied states would be removed from the denominator. 
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Substituting the steady state value of Pi from Eq. ([32]) into Eq. (|1T]) gives: 



1 + 



An 



A2 +4e? 



(44) 



hence we may substitute this into Eq. ()43p to give the final form of the self-consistency condition 
for the Maxwell-Bloch equations: 



(cjo - /X - iK)V'o 



2 , (2ej + iXi_) 

\ ^""^^ WTAiTI(Al7W^ 



(45) 



The laser threshold condition is given by taking ^ — )• in the above equation. If we also take 
9i = 9, = and the usual laser operating condition of Aj_ ^ k one has that lasing occurs at the 
cavity frequency, fi = ujq and the threshold condition has the well-known form: KX±/g'^ = uNq, 
where n is the number of two-level systems. 



Recovering laser limit from non-equilibrium mean-field theory 

This simple laser self-consistency condition can be recovered from equation (|33|) if rather than 
using the frequency dependent forms for F^,b(z^) discussed previously, one instead takes Fa,b 
to be constants . Physically such a limit corresponds to high temperatures. Note that as the 
temperature rises, to keep the bath population fixed, the chemical potential must also vary. We 
will therefore take fj, (x T, and then take the limit T — t- oo. Such a limit has another simple 
interpretation, corresponding to making a fully Markovian approximation, including assuming 
the occupation, as well as the density of states, to be flat, and so writing the Keldysh part of 
the self energy as 



As such, our approach in Eq. (l29l) is Markovian for the density of states of the bath, but 
non-Markovian for the occupation. In terms of quantum statistical (i.e. Heisenberg-Langevin) 
approaches, the distinction is whether the noise should be taken as white noise or coloured noise. 
Assuming the noise correlations to be white, and thus neglecting the frequency dependence of 
occupation, is also the approximation underlying the quantum regression theorem (49|. which 
allows one to relate two-time correlations to the evolution of the density matrix. The role of 
this approximation, and its implications for the fluctuation dissipation theorem are discussed by 



Ford and O'Connel[50| 



If Fa,b are frequency independent, then in Eq. (I33p . the term in the integral proportional to 
1/ will vanish as this is an odd function, and so Eq. ([33|) becomes: 

(wo - M - i«^)^o = Yl diMFB - ^a)^^^^^^- (47) 

Hence, the polarisation of the two-level systems is in this case proportional to the inversion of 
the baths, A'^o = [ub — ua) = —{Fb — Fa)/2 and we have: 



Then, identifying the decay constants in Eq. (1451) as A_l = Ay = 27, Eq. (1481) and Eq. (j45|) are 
equivalent. 
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General properties of mean-field theory away from extremes 

Away from the extremes of laser theory or of thermal equilibrium, the effect of pumping on the 
phase boundary can be understood as a result of competition of two effects: pumping and decay 
add noise, reducing coherence, hence suppressing condensation; on the other hand, for a given 
decay rate, pumping increases the density, favouring condensation. The simplest illustration of 
the first of these is shown in Fig. 10.41 where one sees that as the value of 7 is increased, for a 
fixed K, the critical density required for condensation increases. 




Critical density 



Figure 0.4: Critical temperature as a function of density, showing effects of pumping 
and decay, taking a Gaussian distribution of two-level-system energies with variance 
0.155. Adapted from Ref.^. 



To see the competition between pumping causing dephasing and pumping increasing density, 
one may look at the low temperature limit, shown in Fig. 10.51 plotting the critical value of k as 
a function of 7. Two lines are shown; the solid line has an inverted bath (as would be required 
for the laser limit), the dashed line has a non-inverted bath. In the later case (as illustrated 
in the inset) for small 7, the two-level system energy is too far below the pumping bath, and 
insufficiently broadened by 7, to be populated; for larger 7 the broadening is sufficient, and 
condensation may occur. In the presence of inhomogeneous broadening, the above picture is 
significantly relaxed, since the tail of the density of states can be occupied even if the peak is 
below the chemical potential. 

Low density limit: recovering complex Gross-Pitaevksii equation 

The self-consistency condition of Eq. (j33p can also be related to the idea of the complex Gross- 
Pitaevskii equation providing a mean-field description of a spatially varying condensate. For a 
steady uniform state, the mean-field self-consistency condition may be understood as as (/x^ -|- 
in — u}o)'tpQ = xiV'Oi /^S'lV'Oi where x[V'0)/"5] is a nonlinear complex susceptibility. For a ipo{r,t) 
which varies slowly in space and time [up to an allowed fast time dependence described by a 
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Decay rate, K/g 

Figure 0.5: Critical couplings to pumping bath without inhomogeneous broadening 
and at low temperatures. 



factor exp{—iij,st)], one may consider the local density approximation: 



idt + in 



21 



V 

V{r) - — 
^ ' 2m 



iJo{r,t) = xbPo{r,t)]^po{r,t). 



(49) 



In order to determine the large scale spatial structure of a condensate, or its low energy collective 
modes, it is often sufficient to make a Taylor expansion of the nonlinear complex susceptibility, 
resulting in a complex Gross-Pitaevskii equation: 



— + V{r) + U\iJo\'^ + i [7es{^^B) 
2m ^ 



(50) 



where F represents the simplest form of nonlinearity of the imaginary part, taking a form that 
will ensure stability. 

Depending on the details of pumping included in the model, one may find that by treating 
x[V'(^)] more carefully the susceptibility depends not only on the current value of tp{t), but on its 
history, due to dynamics of the reservoir. [In fact, to correctly reproduce the polariton spectrum, 
one ought to take the excitonic susceptibility to have a resonance at the exciton energy, after 
which a variant of Eq. ()50p . but with the appropriate polariton dispersion will be recovered.] In 
the limit of sufficiently slow dynamics of the system, or when considering steady states, dynamics 
of the reservoir should become unimportant. Results of the complex Gross-Pitaevskii equation 
with or without separate reservoir dynamics may be found elsewhere 4^47] 



5 Fluctuations, and instability of the normal state 



As stated earlier, when introducing the self-consistency condition for the non-equilibrium prob- 
lem, it is not possible to consider minimising free energy when looking at a system coupled to 
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multiple baths, and so it is instead necessary to look for stable steady states. The self consis- 
tency conditions discussed above determine whether a steady state may exist, but not whether 
it is stable. In order to analyse stability, it is necessary to consider fluctuations about a given 
state, and to find whether they grow or decay in time. In addition, the study of fluctuations 
allows one to determine the response functions of the system — in the present context, this 
means the photon Green's function — which in turn will give the physical observables, such as 
photoluminescence and absorption spectra. 

In the non-equilibrium case, both the spectrum of possible excitations (i.e. what is seen 
in the absorption spectrum), and its occupation (i.e. photoluminescence) must be determined 
independently, for which the Keldysh Green's function approach is ideal. In the following, the 
approach to calculating these Green's functions is discussed for both the normal and condensed 
state, and then this approach is applied to understanding the instability of the normal state, 
which allows a clearer interpretation of the relation between the non-equilibrium condensate 
and a simple laser. For the condensed system, the calculations are more complicated due to the 
existence of non-zero anomalous correlations, i.e. {ipk{t)''P-k{t'))', the general structure of the 
spectrum of the non-equilibrium system will be discussed in section [6l 



(51) 



Photon Green's functions in the non-equihbrium model 

To allow for anomalous correlations in the condensed state, it is helpful to write the Green's 
function in a vector space of V'fciV'lfc' J^^t as in the above discussion of the Green's functions 
for the two-level system, this vector space of V'fcjV'Lfc should be combined with the it space due 
to the Keldysh/retarded/advanced structure. Thus, one has four by four matrices, in the basis 

(V'fe,+>^lfc,+,V'fc,_,V'lfc,-)- 

The photon Green's function can be found by solving the Dyson equation, D^^ = Dq^ — S, 
and so to start with, the free photon Green's function is required. The free Hamiltonian in this 
case is just fj^^°^°^ = Yl,k^k'4'\'4^k- In the four by four basis arising from mixing V'fcjV'lfe) some 
elements correspond to Green's functions in which -0, ip'^ are interchanged in order. This means 
that these elements are Hermitian conjugated, giving the form: 

/ u — Cjk — irj \ 

—bj — ojk + iv 

u — bJk + ir] (2i?7)Fo(a; + /i) 

\ -u-Cjk-iT] {2iri)Fo{-u} + ) 

where once again r/ is infinitesimal. In this, we have written all frequencies measured relative to 
//5, meaning that we made the substitution -0^ — )■ e~^^^^^{ipo6kfi + ipk)- 

To this free Green's function one must add self energies arising from two parts of the in- 
teraction Hamiltonian. The first is the coupling between cavity photons and the decay bath; 
the second is the coupling between the photons and the pumped two-level systems. The first 
contribution has a form exactly analogous to the coupling between the two-level systems and 
pumping baths, i.e.: 

p p 
p p 

p p 
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where as before Fij, is the Fourier transform of the 2n^f(ci;) + 1, and the Markovian hmit for the 
bath density of states and couphng constant has been applied to get the final expression; these 
terms thus give a self energy: 



Sdecay('^) 



/ 








+iK 





\ 













— IK 






—iK 





-{2iK)Fi!,{u} + fis) 







V 





+iK 





-{2iK)F^{-uj + fis) 


/ 



(52) 



In calculating the self energy due to the coupling to two- level systems, one may simplify the 
calculation by noting that only S^^^, S^^^^, S^^^, S^^^^ are independent; all other self energies 
can be related to these quantities by Hermitian conjugation and/or swapping uj — )• —lu. To 
generate the diagrams for these self energies, we should first determine the interaction vertices 
that give rise to such self energies. The relevant part of the interaction Hamiltonian here is the 
interaction between two-level systems and incoherent photons, and so the relevant contribution 
to U comes from 



dtH; 



C 



TLS-photon 
int 



dtg 



-oo 
oo 



^{t, f)b]{t, f)d,{t, f) - ^lj{t, b)h\{t, b)d,{t, b) + H.C. 



oo V 2 



'Mt) {h\^{t)d,_{t) + bU€)di^{t) 



+ V-(t)(6l+m+(t) + 6t(t)d,_(t) 



+ H.C. 



(53) 



The self energy diagrams thus consist of diagrams with one incoming and one outgoing photon 
line, connected via the interaction vertices in Eq. (j53p . and the Green's functions for the two- 
level system. As is clear from Eq. (j53p . the vertices all involve the two-level system swapping 
between the a and b states. Just as for the diagrams describing the effects of the bath discussed 
in Sec. O one must also keep track of the it labels on the fields. To calculate, for instance, 
the retarded self energy (i.e. the — -|- component) it is clear that the vertices arising from the 
possible placements of it signs have the form: 

b b 




(any other set of possible it labels on the internal lines will involve a line, and such Green's 

functions vanish). To translate these dia gram s into an equation for the self energy, one must 
use the following Feynman rules (see Refs. 31- 13]): For each interaction vertex there is a factor 
{—ig/^/2), and for each internal Green's function, a factor iG. There is then a prefactor 
where F is the number of closed Fermion loops (F = 1 in the current case), and there is a 
combinatoric factor associated with how the vertices are found from the expansion of U, which 
is the same as in any other diagrammatic approach. Applying these rules, one may write: 



dv 



L "■i'^i ^ 



bib, 



2! \V2j 

For the anomalous case, all that changes is the a, b labels, i.e.: 



z. + a;) + G^ (zv)Gft, 

aj a. ^ ' b' b- 



V + UJ] 



(54) 



(55) 



The component S"*" is just the Hermitian conjugate of S + as above. The component 
vanishes, since it either involves lines, or it involves products of two retarded Green's 
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functions. Since the retarded Green's function is causal — i.e. D^{t, t') oc 6(t — t') — then as a 
function of frequency, all of its poles are in the lower half plane, and so the integral of a product 
of two such functions is equal to zero. ^ The only other surviving component of the self energy 
is thus: 




a a a 



which gives the equation: 



(56) 



As was the case for the retarded components, the only difference between normal and anomalous 
Keldysh components is in the a, b labels, so: 



(57) 



Combining the self energies due to the pumped two-level systems and the self energy due 
to decay with the free inverse Green's function, one can then find expressions for the photon 
Green's functions, and hence observable quantities such as the photoluminescence intensity as 
a function of frequency and momentum, which is given by C{uj) = i (^D^^^ — D^^^ + ^^t |/,) Z^- 
In section [SI the normal state Green's functions are studied: we show how an effective density 
of states and occupation function can be defined, and also show how the behaviour of these 
functions can be related to the structure of the inverse Green's function, and to the stability of 
the normal system. 

In the condensed state, just as in equilibrium, the form of the inverse Green's function 
can be shown to obey the Hugenholtz-Pines relation [58] (see also (sij, Chapter 6]), meaning 
that [D^^^]~^(0, 0) = [D^|^|]~"^(0, 0), which implies there is a gapless spectrum. Just as in 
equilibrium, one may show that the requirement for the Hugenholtz-Pines relation to be sat- 
isfied is equivalent to the mean-field condition, Eq. (|33]l . It is worth noting that as iIjq — )• 0, 
the Hugenholtz-Pines relation (and hence the mean-field condition) become equivalent to the 
condition that: 



P^^t^] ^{'^ = /^eflf, A; = 0) = ^cflf - ^0 + - ^TLsC/^off) = 0, (58) 

for some particular //gfT- (In this expression, the non-condensed self energies have been written 
without the gauge transform of Eq. (I22p . as in the absence of a condensate, there is no reason 
to perform the gauge transformation.) Section [5] will show that the condition in Eq. (j58p also 
corresponds to the point when the normal state ceases to be stable. 



^NB; since the Green's function generically looks like l/o; at large to, the integral of a single retarded Green's 
function depends on the regularisation used. However, for a product of retarded Green's functions, the integral 
is well defined, and so vanishes. 



5. FLUCTUATIONS, AND INSTABILITY OF THE NORMAL STATE 



19 



Normal-state Green's functions and instability 



Focusing on the non-condensed case, the properties of the spectrum are entirely determined by 
three real functions of uj, as one may write: 



D 



R 



D 



-1 



K 



(w) = iC{uj) 



(59) 



The forms of A{oj), B{uj), C{uj) follow from the expressions in the previous section. These some- 
what simplify since we are considering the normal case, and so one has: 



OJ 



^fLs('^) 



D 



-1 



K 



2zkF^(w) - SfLg(tj), 



(60) 



where Srpj^g are the self energies from the pumped two-level systems given by Eq. (iMl) and 
Eq. (j56p . In the following we will first discuss how the forms of A{ijj), B{uj), C{oj) determine the 
spectrum, occupation and stability, and then illustrate this with their forms arising from the 
particular microscopic model discussed above. 

Inverting the matrix of Keldysh Green's functions (using Eq. ([9])), one finds: 



<t^(^) 



A{uj) - iB{uj) 



-iC{uj) 



A(w)2 + 5(a;)2' '^^^^"'^ A(w)2 + 5(a;)2^ 

and then in terms of these quantities, we may write the luminescence spectrum: 

C(w) - 2B{uj) 



2[^(a;)2 + 5(a;)2] 



(61) 



(62) 



Further, by analogy with the equilibrium system, we can explain the form of this expression in 
terms of a spectral weight (density of states) p{oj) = —2'^[D^^ ^{uj)] and an occupation function 



2n^{uj) + 1 = iD^^^{u:)/p{u: 



givmg: 
2B{uj) 



yl(a;)2 + S(L^)2^ 



2B(u) 



1 



(63) 



hence the luminescence is related to these as C{uj) = p{uj)n^{uj) as expected. 

In the absence of coupling to the two-level systems (and hence neglecting Srpj^g in Eq. (f60|l 
, one may clearly identify the role of the three expressions involved here: 

• B{uj) = K is the linewidth of the normal modes 

• A{uj) = uj — (jjk describes the locations of these modes, and 

• Cioj) = 2K(2n^ + 1) describes their occupation. 

However, when coupled to the two-level systems, B{uj) is not a constant, hence firstly, the 
linewidth varies, and more importantly, B{u) may vanish at some value of w. If B(uj) does 
vanish then the occupation diverges, but since the spectral weight vanishes too, the luminescence 
remains finite. 

Physically, this describes the behaviour that would, in equilibrium, be expected at the chem- 
ical potential, as long as the chemical potential lies below the bottom of the band. Note that 
the equilibrium Bose-Einstein distribution diverges at the chemical potential. However, if the 
chemical potential lies below the bottom of the band then the spectral weight is zero at the 
chemical potential and thus the particle number (luminescence) remains finite. 

Out of equilibrium, the system distribution may in general be far from the Bose-Einstein 
distribution. Even so, when near the threshold for condensation, the system distribution shares 
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an important property with the Bose-Einstein distribution: near the frequency where the con- 
densate will emerge [i.e near the point where B(uj) = 0] the system distribution will diverge 
as l/{u} — //eff), just as the Bose-Einstein distribution does. We may thus identify the effect of 
pumping as introducing a chemical potential that has nothing to do with the chemical potential 
of the decay bath. Since B{lj) is given by the inverse retarded Green's function, one may note 
that the inverse Keldysh Green's function does not on its own fix the distribution; it is the ratio 
of Keldysh and imaginary retarded Green's functions that matter. Figure 10.61 shows how the 
spectral weight, occupation and luminescence are related to the zeros of the real and imaginary 
parts of the inverse Green's function. 
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Figure 0.6: Behaviour of inverse Green's functions, and resulting properties of spec- 
tral weight, luminescence and occupation functions in the normal state. Upper panel 
shows the inverse Green's functions (with zeros marked by arrows), and the lower 
panel shows the various physical correlations of interest. Adapted from Ref . (29| . 



Zeros of A{uj),B{uj) and stability 



Although a zero of B{uj) alone does not cause the luminescence to diverge, a simultaneous zero 
of A{uj) and B{u}) will. The stability of the system can be seen to change when this occurs, as 
will be discussed next. When near a simultaneous zero, one may expand A{uj) = a{uj — and 
B{uj) = j3{oj — /ieflf), and so: 



I , [ ("C + «/?/^efr)(a 
(a -I- ip) oj 



a2 + /?2 



hence the actual poles are at frequencies: 







a2 +/32 



(64) 



(65) 



These poles determine the time dependence of the retarded Green's function, so if > then 

the pole has the wrong sign of imaginary part and the normal state is unstable. When yUgg = ^, 

-1 



then this means there is a value u 



^ for which 



[uj^k = 0) = 0, which as 
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discussed in Eq. ()58p is equivalent to saying that the mean-field consistency condition can be 
satisfied. Hence, instability of the normal state, and the existence of a condensed solution will 
occur together. 

It is helpful here to explicitly write B{uj), in order to understand the origin of its zeros, and 
what parameters determine their location. In the non-condensed case, the fermionic Green's 
functions that come from inverting Eq. (|30l) have a simple form: 



and so substituting these into Eq. ([SI]) , and taking the imaginary part one may write: 

J ^'i [(i' + a, - £i)2 + 72J [(i' + eif+ 72 

For -6(0;) to have zeros, it is necessary that the second term (which originates from pumping) 
should be negative, and should overcome the first term (which originates from decay). With 
Fa,b{^) = tanh [li{v it /iB/2)/2], it is clear that this criterion requires to be sufficiently large. 
As such, the following scenario describes what happens as ^ib is increased: 

Very weak pumping. For large negative ^b^ one finds that Fb{v + o;) — -Fa(i^) is always 
positive, and so no zero of B[ijj) exists. 

Subcritical pumping. For less negative values oi ^b, there is a range of uj for which B{uj) is 
negative, indicating a range of gain in the spectrum. The boundary of this region, where 
B{uj) = defines an effective chemical potential ^eff; but since /^eff < C the normal state 
remains stable. 

Critical pumping. At some value of /is, one finds that /Xeff = ^, meaning that at this value of 
= f^cfi = one has D^^_^{uj*) = 0. Hence, the gap equation first has a solution at this 
point, there is a real divergence of the luminescence, and the normal state is marginally 
stable. 

Supercritical pumping. Above this critical value oi ^b-, the normal state would have ji^s > 
and so would be unstable. 

The actual behaviour for the polariton model of Eq. ([1]) is shown in Fig. 10. 7t one can see that 
a pair of zeros of the imaginary part emerge, and then one crosses the bottom of the polariton 
modes. Note that in equilibrium, we have fieS = /i = /^b at all conditions, and so only the 
last three stages of the above scenario exist; condensation occurs when the chemical potential 
reaches the bottom of the band. It is also important to note that the above scenario means that 
the Bose-Einstein distribution is not the only distribution that would allow condensation. Any 
distribution which has the above property, i.e. a divergence at some frequency for given values of 
the control parameters (density, coupling constant, etc.), is sufficient for quantum condensation 
in bosonic systems. 



Simplified form of distribution function in high-temperature limit 

The way in which the effective distribution is set by the balance of pumping and decay can be 
demonstrated more clearly by specialising to the case of 7 ^ T, for which the pumping bath 
occupation functions do not change significantly across each Lorentzian broadened peak (but 
may vary between the two peaks). In addition, consider taking gi = g,ei = e, so that sums of 
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Figure 0.7: Variation of energies of zeros of real part of inverse Green's function 
and imaginary part as density is varied via chemical potential of pumping bath. 
The three solid lines correspond to (from the bottom) lower polariton, exciton, and 
upper polariton respectively. The point where the dashed line crosses solid line is 
where the condensation occurs. Adapted from Ref [l^. 



two-level systems can be replaced by factors n. Then the expression in Eq. ()67p can be simplified 
to give: 

B{uj) = K + ng'j '^^'^' 7\ (68) 



^^ jFB{e)-FA{e-u )] 



where we have performed the integrals assuming the distributions are effectively constant. ^ By 



applying the same approach to [D one has: 



2 

K{2n^{u;) + 1) + ^^2^xt[1 " ^B(e)i^A(e - u:)] 

2n,iu) + 1 = iu.-2er + i,- _ 

+ (.-2"y+47^ [^"(^)-^^(^-")] 

^From this expression one may first note that if 7 = (or more generally if k ^ g'^j/[{LU — 2e)^ + 
47^]), the the system distribution is the same as the distribution of the decay bath (the photons 
outside the cavity) and so n^{uj) = n^(a;). On the other hand, if k = 0, (or more generally, if 
K <C g^^/[{uj — 2e)^ + 47^], which can occur near a; = 2e), the distribution is set by the pumping 
bath. In this case, the important terms in Eq. ([69|) are: 



2n^{uj) + 1 = ^_,X T^i = coth e-^-e + w- ^ ), (70) 



FB{e)-FA{t-uj) \2 



which is a Bose distribution with the temperature and chemical potential of the pumping bath. 
Thus, the photon distribution interpolates between the decay and pumping bath, depending on 
the efficiency of coupling as a function of energy. An illustration of how this might look when 
the chemical potential of the decay bath is not too dissimilar from the pumping bath is shown 
in Fig. 10.81 however for realistic parameters, the chemical potential of the decay bath should be 
taken to —00. 



^ Formally, the approximation consists of performing the contour integral, taking into account the poles at 
1/ = —iu + ei+i^ and u — —£ + 17, but neglecting the poles from Fa. nil') which are ai v — {—lo + ij.b/2, — /is/2} + 
j(2n + l)nT, along with neglecting in evaluating the residues. 
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Figure 0.8: Cartoon of occupation function set by competition of bosonic bath 
and fermionic bath, with effect of fermionic bath moderated by a Lorentzian filter 
depending on excitonic energy. The chemical potential of the decay bath is at —9, 
and that of the pumping bath is just below zero. 



Normal-state instability for a simple laser 

As for the mean-field theory, it is instructive to compare the results of Sec. \E\ to those for a 
simple laser, in which pumping tries to fix the inversion of the gain medium, independent of 
frequency. The instability of the normal state can still be determined by the inverse retarded 
Green's function, which may in turn be found by the response of Eq. (|39l )- ()4ip to an applied 
force _Pe~*'^* acting on the photons. If the force is weak, then Eq. (|^T|) reduces to Ni = Nq, and 
taking A_l = 27 as found previously, the equations to solve are: 

dtip = -iujQip - + X] 9iPi + Fe-'''\ dtPi = -2ieiPi - 2jPi + gi;No, (71) 

i 

hence writing the response as tp = iD^^^{uj)Fe~^'^^ , and eliminating Pi gives: 

[««,i-M=.-„„+,.+x:-^-|^. (72) 

i 

As in the mean-field case, this same equation can be recovered from the microscopic non- 
equilibrium model by taking Fa^b to be independent of frequency, and identifying A^o = —{Fb — 
Fa)I2. The form of the inverse retarded Green's function makes much clearer the implications 
of this absence of frequency dependence. For the imaginary part of Eq. ()72p to be zero, it is 
clearly necessary that > 0, so a region of gain can only exist when inverted. 

In the special case of = e = ujQ/2,gi = g, the zeros of the real and imaginary parts can be 
found explicitly to be 



^l,s = 2e±^g^nNo^-4J''., ^ = 2e, 2e ± ^"47^ " 9^Non, (73) 

where n is the number of two-level systems as before. From the zeros of the imaginary part, 
one sees that a region of gain exists only for Nq > 2^7/(7^ n (note that this is the laser threshold 
condition discussed in section U]). On the other hand, a splitting of the zeros of the real part 
^ exists only if Nq < —Aj'^/g'^n. Thus the instability of the normal state only occurs after the 
normal mode splitting has collapsed. This is illustrated in Fig. 10.91 In this figure, it is also 
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clear that as soon as there is a region of gain, there is an instability. This is quite different 
from Fig. 10.71 where a region of gain, and thus zeros of the imaginary part, emerged at a lower 
pumping strength than was required for the instability. This meant that in the non-equilibrium 
condensate, a diverging distribution function exists before condensation occurs, whereas for 
Fig. 10.91 the distribution function has no divergence in the normal state. ^ 
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Figure 0.9: As for Fig. 10.71 but for the results of the Maxwell-Bloch equations, 
showing the rather different behaviour in the extreme laser limit, plotted for wq = 2e. 



6 Fluctuations of the condensed system 

When condensed, the derivation of the spectrum from the inverse Green's function written 
previously becomes much more involved, but the essential features of the spectrum can be 
determined by considering the symmetries that the system must possess — the results of this 
analysis are confirmed by the exact expressions for the inverse Green's functions. In particular, 
one may combine the Hugenholtz-Pines relation, mentioned at the end of Sec. El with the 
analytic properties of the Green's functions which imply [Z)^^|]~-^(a;,p) = [-D^|^]~-^(— a;,p)*, 

= ^From these general considerations, one may find that the 

most general structure for sufficiently small a;, /c is: 

C C 

D^t^^, k) = = w2 + 2ia;x-c2A;2' ^^^^ 

where x is an effective linewidth, and c an effective sound velocity. The form of this expression 
is dictated by: the need to combine symmetry under k — )• —k; the existence of a finite linewidth; 
and the pole at cj = 0, A; = that is ensured by the Hugenholtz-Pines relation. Higher order 
contributions could exist (and in fact do exist) for larger a;, k, but the w. A; — )■ structure is fixed 
by these considerations. 

The above structure means that the poles of the Green's function for small k are diffusive, 
i.e. OJ* = —ix±iVx'^ — c^k"^, meaning that long wavelength excitations decay, but with a lifetime 



''if one considers the more general case with detuning, ujo 7^ 2e, a region of gain may appear before the 
instability occurs. Furthermore, if one also has inhomogeneous broadening, e 7^ ej, and different inversion for 
different two-level systems, a region of gain can coexist with a splitting of the normal states. Ifowever, the results 
for the non-equilibrium condensate shown in Fig. 10. 71 had neither detuning nor inhomogeneous broadening; hence 
in the absence of such complications, the difference between the non-equilibrium condensate and a simple laser 
are particularly obvious. 
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that diverges for one mode as A; — t- 0. This same form is also recovered from other approaches to 



non-equilibrium condensates |45l|. including also the case of a parametrically pumped polariton 
system 60(] . Note that if one were to naively extract a Landau critical velocity from the real part 
of (jj* , then this critical velocity would vanish. There has been some work on how the concept of 
the Landau critical velocity may be generalised for parametrically pumped condensates [6 ll-t63l] . 
however the full implications of the diffusive structure on superfluidity of incoherently pumped 
non-equilibrium condensates remains an open question. Because the polariton system is two- 
dimensional, phase fluctuations can be expected to play a particularly important role, therefore 
the remainder of this section will discuss how the above form of the Green's function determines 
the long-time correlations, and hence the lineshape, and how this connects to other approaches 
to deriving the polariton lineshape. 

To take full account of the phase fluctuations, one must reparameterise the fluctuations as 
ijj = \/ p + TTe^'^. In order that one works with fields for which there is a macroscopic expectation 
of this reparameterisation must be performed in real space, and in terms of the fields on the 
forward and backward contours, rather than the symmetric and antisymmetric combinations 
of these fields. (Note that the macroscopic expectation of the anti-symmetric combination 
vanishes 0].) To describe the long-time correlations, we wish to find the first order coherence 
function D^^^{t) = — i(Tc['0(t, /)'0^(O, 6)]), and corresponds to the Fourier transform of the 
luminescence spectrum, C{u}). Since it is the phase fluctuations that dominate the long time 
behaviour, one may write this asymptotic behaviour in the form: 

Kh^^^ ^ P<^^ ^^^P ('^(^^ - = -"0^ exp[-/(t)], (75) 

where pqc is the quasi-condensate density. The function f(t) is given by the phase-phase 
correlation functions, and in two dimensions is given by: 

m = ^ [dZH) - DiliO)] =l^l^[l- e-'] ^I?£(-, k). (76) 

Note that expressions (j75|) and ([76l) are determined by taking the phase fluctuations to all orders. 
The density fluctuations give no time dependence at long times, their effect appears only in the 
difference between the quasi-condensate density pgc and the total density p. 

Since corresponds to the luminescence spectrum, its relation to Keldysh and retarded 
Green's functions is as in Eq. (|62|) . Assuming that the condensation arises due to pumping, 
then as in Sec. [5l the frequency dependence near the effective chemical potential arises from the 
behaviour of the inverse retarded Green's function — the frequency dependence of the inverse 
Keldysh Green's function has no particular singularities near this point. In this case (which is 
also what is found from the full calculations of the microscopic theory) , the singular behaviour 
of the L>(^ is given by oj^^ ~ l-Df, and so: 



duj_ r kdk jcyp) [1 - e- '^*1 

2tt J 2tt luj"^ -\-2iujx - c^k'^l 



m =1^1 ^^.r^ (77) 



(The factor of 1/p occurs from the relation of phase-phase Green's functions to V'lV'^ Green's 
functions). As one expects for a two-dimensional system, after integrating over uj, the above 
integral reduces to an expression ~ J dk/k, and so one has logarithmic behaviour, cut off at high 
A; by a maximum energy of excited modes, and at small k by the time dependence. At small k, 
the poles of the ui integral are at w = ±2ix,±i{ck)'^ /2x; the flrst of these has a finite residue 
as — )• 0, while the latter has a residue that is is diverging, and thus dominates the behaviour. 
The asymptotic behaviour is thus given by: 

/(t) = [ ^ [l - e--''='*/2xl (78) 

■'^ ' J 2-11 4x(cA;)2 L J ^ ^ 
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(where C is a new constant). This expression has a cutoff for small k given by /c ~ \Jx/t/c; 
thus one still has power law correlations as for an equilibrium two dimensional gas, but with a 
different power law, now set not only by the condensate density but also by the pumping and 
decay strength. Were one to calculate also the long-distance correlations at equal times, one 
would note another difference from equilibrium. In equilibrium, the decay of long-time equal- 
position correlations, and long-distance equal-time correlations have the same power laws. For 
the spectrum in Eq. ([7H) . the power-law for long-distance equal-time correlations is twice that 
of long-time equal-position decay. This is because the low momentum cutoff for long distances 
is always A; ~ 1/r, whereas the long-time cutoff is A; ~ 1/ci in equilibrium, but k ~ \/x/t here. 

Finite-size effects — lineshape of trapped system 

For a confined system, the integral over k modes is replaced by a sum over a discrete set of 
modes; i.e.: 

In this form, one may then consider how the value of the sum depends on the relative size of the 
mode spacing AE', the low energy cutoff s/x/t^ and the maximum energy i^max- Let us assume 
the maximum energy is large, then we have a picture something like Fig. 10.101 
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Figure 0.10: Spacing of discrete energy levels, and upper/lower cutoff energies 



The sum can be divided into parts above and below the low energy cutoff, giving: 

1 



f{t) ^ C 



n=0 



t 

8^ 



(80) 



If both of these sums have many terms, then they may be approximated by integrals, and if 
the density of states z/(^) is uq£^ as it would be for a two-dimensional system with ^„ = cpn (as 
illustrated in Fig. lO.lOp then this becomes: 



f{t) ^ CVo 



/xjt 



8x 



4xe 



/x/i 



16x 4x I V X 




(81) 



What is to be noted here is that the number of terms in the first part compensates the t 
dependence, leading to a harmless constant. If however the number of terms in the first term is 
small, or is in fact truncated at its minimum value of one (which will inevitably occur for large 
enough t), then one instead has something of the form: 
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and so one has exponential decay of coherence at long times, arising from the restricted number 
of modes. 

If the mode spacing is such that only a single mode is involved, then Eq. (j79p reduces to 
a single term in the sum, = 0, and the result becomes very similar to the form found from 
models of phase noise for a single mode condensate [64, 65], for which the lineshape interpolates 
between Gaussian and Lorentzian: 

hence the decay of coherence varies from at short time (giving Gaussian lineshape at high 
frequencies) to t at long times (giving a Lorentzian peak at low frequencies). This result is 



exactly as one expects for phase noise from varying densities 65|: 



dttj) = -iUN, dtN = -TN + Fr(t), (84) 

where Fr is a Gaussian delta correlated noise noise noise with strength Pp- Solving these 
equations in Fourier space, one has: 

».P> = 5(|iV„ft = ;;,(5^, (85) 

which is the same form as in Eq. ()83p 

This section thus shows another distinction between condensates and lasers in terms of many- 
mode or single mode fluctuations. If the system is large the spatial fluctuations resulting from 
the continuum of modes give rise to (in two dimensions) a power-law decay of correlations as 
for an infinite, equilibrium, two-dimensional quasi-condensate. For smaller systems, or at longer 
times, the power law crosses over to exponential decay, given by a fluctuations within the single 
lowest energy mode (the other modes are too high in energy to be relevant), as is characteristic 
for lasers. 



7 Summary 

This chapter has discussed in detail how the non-equilibrium Green's function formalism can be 
applied to study a model of microcavity polaritons, driven out of equilibrium by coupling to two 
baths. This model system, while not incorporating all features of the real system, allows one to 
make particularly transparent connections between laser theory and equilibrium descriptions, as 
well as allowing clear illustrations of the consequences of the approximations typically used for 
simple lasers. By considering steady states of the system in which there is a coherent photon 
field, one finds a criterion for condensation to occur, and can find a self-consistency condition 
which determines how the amplitude and frequency (effective chemical potential) of the coherent 
field depend on the strength of the pumping and decay. By considering fluctuations about steady 
states, one can determine whether a given steady state is stable, find the spectrum of possible 
excitations, and find how this spectrum is populated. 

Starting from the normal state, without a condensate, and increasing pumping strength, 
one finds that fluctuations about the normal state become unstable at the same point that a 
condensed solution appears. The scenario by which this instability occurs on increasing pumping 
strength is quite instructive. As pumping strength increases, a region of energies for which 
there is gain appears in the spectrum. The energy dividing this region of gain from regions of 
loss defines an effective chemical potential, at which the non-equilibrium distribution function 
diverges. Instability occurs at a higher pumping strength, when this effective chemical potential 
(and thus the region of net gain) reach the normal modes of the strongly coupled system, at 
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which point polariton condensation occurs. Such a description unites the lasing picture of gain 
exceeding loss with the equilibrium picture of the chemical potential reaching the bottom of the 
band. 

While the above description allows polariton condensation to be discussed in the language 
of laser theory, the results are rather different from the normal limits assumed for a simple 
laser theory. However, simpler laser theory results can be recovered within the model discussed 
here, as corresponding to a high temperature limit. In this high temperature limit, pumping 
corresponds to effectively white noise, and this was shown to mean that gain only exists when 
pumping bath is inverted. This has the consequence that in this high temperature limit, las- 
ing and strong coupling do not coexist, whereas they can in the low temperature polariton 
condensate. 

When considering fluctuations about the condensed state, a somewhat different distinction 
between simple lasers and the polariton condensate emerges: the effect of finite system size, and 
the spectrum of collective phase modes. For an infinite two dimensional system, the decay of 
coherence at long distances and long times is power law, as in equilibrium (but with different 
powers). For a finite system, the effects of finite lifetime and finite size combine to lead to 
exponential decay at long times. In the limits of very small system size, the standard result for 
phase noise in a single mode condensate is naturally recovered. 

To summarise, the approach presented here provides a way to connect a number of different 
approaches to equilibrium and non-equilibrium condensates, as well as theories of lasers, in a 
transparent manner, allowing one to understand the significance of various approximations, as 
well as the relations between some of the other approaches one may use. 
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